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Prefacio 


Hace unos afios, la Universidad Social me invitó a dictar 
dos conferencias sobre el método de Montecarlo para los 
oyentes de la Facultad de Técnica de cálculo. Yo acepté, 
Poco antes de la primera conferencia me onteré con espanto 
de que los oyentes, en su mayor parte, ignoraban la Teoría 
de las probabilidades. No había modo de retrnceder. Sobre 
la marcha hubo que incluir en la conferencia un apartado 
adicional con el fin de exponer a los oyentos las nociones 
principales de la Teoría de las probabilidades. 

БІ contenido de estas conferencias, dictadas а lo largo 
de varios afios, con el tiempo iba estabilizándose. En la 
edición presente se ha conservado también el apartado adi- 
cional (§ 2) al que quiero dedicar unas palabras. 

Cualquier persona habrá empleado alguna vez las pala- 
bras «probabilidad» y «variable aleatoria». La idea intuitiva 
de la probabilidad (en tanto que frecuencia) corresponde 
más o menos al contenido verdadero de este concepto. En 
cambio, la idea intuitiva de la variable aleatoria está muy 
lejos, como regla, de Ja definición matemática. Por eso, en 
el $ 2 damos por conocido el concepto de probabilidad y 
explicamos sólo ol concepto más complejo de variable alea- 
toria. Este parágrafo no puede reemplazar el curso de la 
Tooría de las probabilidades: la exposición es simplificada 
y no trae demostraciones. Sin embargo, permite hacerse 
una idea sobre las variables aleatorias suficiente para 
comprender los aspectos elementales del método de Monte- 
carlo. 

El objetivo principal de este libro es sugerir a los espe- 
cialistas de las más diversas ramas la idea de que en el 
campo de sus actividados existen problemas que se pueden 
resolver por el método de Montecarlo, 

Los problemas que se analizan en las conferencias son 
variados y suficientemente simples. Sin embargo, no pueden 
abarcar, claro está, todas las esferas de aplicación del mé- 


todo. Me limitaré a dar un ejemplo. En el libró no se habla 
para nada de la medicina. No obstante, los métodos del $ 7 
permiten evaluar las dosis en la radioterapía. Teniendo un 
programa de cálculo de los rayos absorbidos por distintos 
tejidos del cuerpo, se puede dosificar y orientar la radiación 
del modo más eficiente cuidando de no dafiar los tejidos 
Banos, 

Este libro abarca todo cuanto se trató en las conferen- 
cias; contiene una exposición más detallada de los ejemplos; 
además, se han agregado los parágrafos 9 y 10. 


1. Sóbol 


INTRODUCCION 


[zn 
GENERALIDADES 
DEL MÉTODO 


E! método de Montecarlo es un método numérico que 
ermite resolver problemas matemáticos medianto la simu- 
ación de variables aleatorias, 

1.1. Orígenes del método de Montecarlo. Se considera 
como fecha de nacimiento del método de Montecarlo el año 
1949 en el que apareció el artículo titulado «The Monte Carlo 
method»). La creación de este método suele ligarse a los 
nombres de los matemáticos norteamericanos J. von Neu- 
mann y S. Ulam. En la Unión Soviética los primeros artícu- 
log dedicados al método de Montercarlo aparecieron en 
1955 y 1956 3). 

Es curioso que la base teórica del método era bien conoci- 
da desde hace mucho tiempo. Es más, algunos problemas de 
la Estadística se resolvían a veces empleando las muestras 
aleatorias, o sea, aplicando de hecho el método de Montecarlo. 
Sin embargo, hasta la aparición de las máquinas calculadoras 
electrónicas (MCE), este método no encontraba aplicaciones 
suficientemente amplias ya que la simulación a mano de 
variables aleatorias constituye un proceso muy laborioso. 
Es decir, la aparición del método de Montecarlo en tanto 
que un método numérico de gran universalidad so hizo po- 
sible sólo gracias a la creación de las MCE. 

El nombre de «Montecarlo» se debe al de una población 
del principado de Mónaco, celebre por su casa de juego. 
Resulta que uno de los aparatos mecánicos más sencillos 
que permite obtener variables aleatorias es ... la ruleta. 
Acerca de esto hablaremos en el $ 3, En cambio, conviene, 
por lo visto, responder aquí a una pregunta frecuente: 
«¿Ayuda o no el método de Montecarlo a ganar еп el juego 
de la ruleta?» No, no ayuda. E incluso no tiene nada que 
ver con este juego. 

1,2. Ejemplo. Para que el lector se haga una idea más 
clara de las cuestiones que vamos a tratar, consideremos 
un ejemplo muy sencillo. Supongamos que debemos calcular 

3) Metropolis N., Ulam S., The Monte Carlo me- 
thod, J. Amer. statistical assoc., 1949, 44, № 247, 335—341. 

2) Se trata de los artículos de У, V. Chavchanidze, Yu, A. Shréi- 

der y V. S, Vladímirov. 


2-216 


10 


el área de una figura plana S, Esta puede ser una figura 
arbitraria con una frontera curvilinea, definida grafica o 
analiticamente y compuesta de uno o varios pedazos. Supon- 
gamos que se trata de la figura representada en la fig. 1 


ка 


y supongamos que toda la figura está comprendida dentro 
de wn cuadrado de dimensión 1. 

'Tomemos en el cuadrado un nümero ;V de puntos alea- 
torios. Sea М” el número de puntos que aparecen dentro de S. 
Es obvio, por razones geométricas, que el área de 5 es igual 
aproximadamente а la razón N'/N. Cuanto mayor sea № 
tanto mayor será la exactitud de esta estimación. 

En el ejemplo representado en la fig. 1 se ha escogido un 
total de № = 40 puntos. De éstos, М” = 12 puntos aparecen 
dentro de 8. Tenemos N'/N = 12/40 = 0,30 mientras que el 
valor exacto del área de S es 0,35 !). 


1) En la práctica el método de Montecario no se 
aplica al cálculo de ároa de figuras planas ya que existen para ello 
otros métodos que, a pesar de ser més complejos, garantizan una 
exactitud mucho mayor, (Siguo en la pág. 11). 
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1.3. Dos peculiaridades del método de Montecarlo. 
La primera consiste en que su algoritmo tiene una estrue- 
tura muy sencilla. Como regla, se elabora primero un pro- 
grama para la realización de una prueba alestoria (en el 
ejemplo del punto 1.2 hay que escoger un punto aleatorio 
en el cuadrado y comprobar si pertenece o no a $). Después, 
esta prueba se repite N veces de modo que cada experimento 
sea independiente de los restantes y se toma la media de 
los resultados de todos los experimentos. 

Por esto el método de Montecarlo se denomina a veces 
método de pruebas estadísticas. 

La segunda peculiaridad consiste en que el error es, como 
regla, proporcional a la magnitud V D/N, donde D es una 
constante y N es el número de pruebas. Esta fórmula per- 
mite ver que para disminuir el error en 10 veces (en otras 
palabras, para obtener en el resultado otra cifra decimal 
exacta) es preciso aumentar N (o sea, el volumen del trabajo) 
en 100 veces. 

Se pone de manifiesto la imposibilidad de alcanzar en este 
camino una elevada exactitud, Por eso suele decirse que el 
método de Montecarlo resulta especialmente eficaz en la 
solución de problemas en les cuales se necesita conocer el 
resultado con poca exactitud (del 5 al 10 por ciento). 

Sin embargo, un mismo problema puede ser resuelto 
aplicando distintas variantes del método de Montecarlo 5) 
a las que corresponden diferentes valores de D. En numero- 
sos problemas se logra elevar considerablemente la exacti- 
tud escogiendo un procedimiento de cálculo al que le co- 
rresponde un valor mucho menor de D. 


1.4. Problemas que permite resolver el método de Mon- 
teearlo, En primer lugar, el método de Montecarlo permite 
simular cualquier proceso cuya marcba depende de factores 


Sin embargo, el método de Montecarlo expuesto en nuestro ejem- 
plo permite con la misma facilidad calcular el «volumen multidimen- 
sionals de un cuerpo en un espacio multidimensional. Es más, sucede 
a menudo que el método de Montecarlo es en este caso el único método 
numérico que permite resolver ol problema. 

2) En la literatura especial con frecuencia se suelo hablar aho- 
ra de los métodos de Montecarlo (en plural) para subrayar con 
ello que un mismo problema se puede resolver mediante la simulación 
de distintas variables aleatorias. 


ГА 
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aleatorios. En segundo lugar, en muchos problemas mate- 
máticos, que no tienen la menor relación eon cuestiones alea- 
torias, se puede inventar un modelo probabilístico artificial 
(e incluso más de un modelo) que permite resolver estos 
problemas. En realidad, esto es lo que hemos hecho en el 
ejemplo del punto 1.2. 

Por consiguiente, se puede hablar del método de Mon- 
tecarlo como de un método universal en la solución de proble- 
mas matemáticos. 

Es de especíal interés el hecho de que en algunos casos 
conviene abandonar la simulación del proceso aleatorio real 


Fig. 2 


y concentrarse en el análisis del modelo artificial. Una situa- 
ción de este tipo es considerada en el $ 7. 

1.5. Algo más sobre el ejemplo. Volvamos al ejemplo del 
punto 1.2. Para poder realizar el cálculo debemos escoger en 
el cuadrado puntos aleatorios, ¿Cómo realizar esto? 

Imaginemos el siguiente axperimento. La fig. 1 (aumen- 
tada proporcionalmente), con la figura S y el cuadrado, se 
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coloca sobre un muro y sirve de blanco para un tirador que, 
situado a cierta distancia, dispara JV veces apuntando al 
centro del cuadrado. Claro está que no todas las bales darán 
exactamente en el centro: perforarán el blanco en № puntos 
aleatorios 3). ¿Es posible estimar el área de 5 a partir de 
éstos puntos? 

En la fig. 2 viene representado el resultado de un experi- 
mento semejante. En este experimento se tiene N =40, 
N' = 24 y N'IN =0,60 lo que representa casi el doble del 
valor exacto (0,35) del área. Por supuesto, tratándose de un 
tirador de gran puntería los resultados del experimento serán 
muy malos ya que casi todas las balas pegarán cerca del 
centro y perforarán 57). 

Es fácil persuadirse de que el método que hemos aplicado 
para calcular el área tendrá validez sólo si los puntos alea- 
torios, además de ser aleatorios, estarán «distribuidos uni- 
formemente» sobre el cuadrado. Para dar un sentido exacto 
a estas palabras es preciso estudiar la definición de variable 
aleatoria y algunas de sus propiedades. Estos conocimientos 
se exponen, sin demostración, en el $ 2. El lector que haya 
leído cualquier curso de la Teoría de las probabilidades puede 
omitir por completo el $ 2 a excepción del punto 2,5. 


>) Aceptamos que el tirador está lejos de ser el 
campeón del mundo y se encuentra в una distancia suficientemente 
grande del blanco. 
2 En el punto 4.5 explicaremos cómo han sido escogidos los 
puntos aleatorios en las figs. 1 y 2. 


CAPITULO 1 


SIMULACION 
DE VARIABLES ALEATORIAS 


$2. 
VARIABLES ALFATORIAS 


Suponemos que el lector está más o menos familiarizado 
con el concepto de probabilidad y pasamos directamente al 
copcepto de variable aleatoria. 

Las palabras «variable aleatoria» suelen emplearse para 
subrayar que se ignora el valor concreto que tomará esta 
variable. Sucede también que estas palabras encierran inclu- 
so un desconocimiento de cuál es esta variable, 

En cambio, el matemático emplea estas mismas pala- 
bras «variable aleatoria» dándoles un contenido positivo 
bien determinado. Efectivamente, dice él, no conocemos 
el valor que tomará esta variable en un caso concrelo dado 
pero sabemos qué valores puede tomar y también las proba- 
bilidades de unos u otros valores, A base de estos datos no 
podemos predecir con exactitud el resultado de una prueba 
relacionada con esta varisble aleatoria pero sí podemos 
prever con gran seguridad los resultados de un gran número 
de pruebas. Cuanto mayor sea el número de pruebas, mayor 
exactitud tendrán nuestras predicciones. 

Es decir, para definir una variable aleatoria es preciso 
señalar los valores que puede tomar y las probabilidades 
de estos valores. 

2.1. Variables aleatorias discretas. Se dice que la variable 
aleatoria Ё es discreta si puedo tomar un conjunto discreto 
de valores z,, 23, .. лї}. 

Una variable вазона, "discreta E se define mediante la 


tabla 
Т) Ly eee Tp ) 
= T 
M Paves Pa 1" m 
donde x, 2), . . ., z4 son los valores posibles de la variable 
ЕУ Рі» Pa ++ +» Pn Son las probabilidades que les corres- 


ponden. Hablando con más precisión, la probabilidad de 


1) La Teoría de las probabilidades ostudia también 
variables aleatorias discretas que pueden tomar un número infinito 
de valores. 
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que la variable aleatoria Е tome el valor т; (designaremos 
esta probabilidad por P (E = z,)) es igual a pr: 


P (6 = а) = pr 


La tabla (T) se denomina distribución de la variable alea- 
toria, 

Hablando en términos generales, los números z,, Za, . + 
«as т, pueden ser cualesquiera. En cambio, las Probabili- 
dades ру, Pas . . ·, Pn deben cumplir dos condiciones: 

a) todos los números p; deben ser positivos: 


рі > 0; (1) 
b) la suma de todos los p, debe ser igual a 1: 
Pictpac eb рь = 1. (2) 


La última condición significa que E debe зай 
tomar en cada caso uno de los valores ғу, 2), .. 

Se denomina esperanza matemática de la variable illoa toris 
£ el número. 


ME= 5 р. (3) 
із 


Para aclarar el significado real de esta magnitud, la 
representaremos en la forma 


м-2 np 22 Pr 


De aquí puede verse que МЕ es el valor medio de la varia- 
ble E con la particularidad de que los valores z; de mayor 
probabilidad figuran өп la suma con pesos mayores 1). 


1) La media ponderada aparece con frecuencia en 
las más diversas ramas de la ciencia. Por ejemplo, en la Mecánica: 
si en los puntos а, Lay +++ + Zn (del eje Oz) están concentradas las ma- 
SAS ті, May eee "laa abscisa del centro de gravedad de este sistema. 
se determina mediante la fórmula 


3 (Ў Mj. 


Por supuesto, la suma de todas las masas puede ser en esto caso distinta 
de la unidad. 
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Señalemos las propiedades principales de la esperanza 
matemática: si c es una variable no aleatoria, se tiene 


M(E +e) = М + с, (4) 

M (ct) =cME; (5) 

si E y т son dos variables aleatorias cualesquiera, se tiene 

M (E + n) = ME + Mn. (6) 

Se denomina varianza de la variable aleatoría & el nü- 
mero 

DE = М (Е — МЕЗ. (7) 


Por consiguiente, la varianza DE es la esperanza matemáti- 
ca del cuadrado de la desviación de la variable aleatoria E 
de su valor medio МЕ. Es obvio que siempre se tiene DE > 0. 

La esperanza matemática y la varianza son las caracterís- 
ticas numéricas más importantes de la variable aleatoria E. 
¿Qué valor práctico tienen? 

Supongamos que hemos observado varias veces la va- 
triable E obteniendo los valores E,, En, - + ., Ey (cada uno 
de los cuales es igual a uno de los números тү, 2а, - - ., zz 
entonces la media aritmética de estos valores será próxima 


a ME 
Ath +e. Ы) л ME 9) 


y la varianza ОЕ caracterizará el esparcimiento de estos 
valores alrededor del valor medio ME. 
Empleando (4), (5) y (6) podemos dar otra forma a (7): 


Di = М Iz? — 2М&-&-+ МӘЗ = M (59) — 
— 2M5-M5 + M$)’, 


05 = М (E?) — (MB 9) 


Resulta más fácil, como regla, calcular la varianza apli- 
cando la fórmula (9) en lugar de la fórmula (7). 

Señalemos las propiedades principales de la varianza: 
Si e es una variable no aleatoria, se tiene 


D (E + c) = DE (10) 
D (c£) = c*DE. ш) 


de donde 
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En la Teoría de las probabilidades desempeña un papel 
importante el concepto de independencia de variables alea- 
torias. En realidad se trata de un concepto bastante com- 
plejo; pero en los casos más sencillos Ja cuestión se aclara 
fácilmente !). 

Para dos variables aleatorias independientes E y т son 
válidas las relaciones 

M (En) = ME-Mn, (12) 
D (È + n) = DE + Dn. (13) 

Ejemplo. Consideremos la variable aleatoria x con la distribu- 

clón siguiente 
1 2 3 4 5 8 

к= (us 1/6 1/6 1/6 4/6 "E 

Es obvio que podemos realizar esta variable considerando el nüruero de 

puntos que aparece al arrojar un dado: todos los valores son igual- 


mente probables, Caloulemos la esperanza matemática y la varianza 
de x. Según la fórmula (3), tenemos 


+ 
Muze dep 65 


1 


=3,5 


Según la fórmula (8), tenemos 
Ох= M (2) — (Mao 2 э ва. (ам 2,047. 


Ejemplo, Consideremos la variable aleatoria 0 con la distribución 


Siguiente 34 
o=(y, 4): 


Podemos realizar esta variable considerando el juego de cara y cruz en 
el quo a la cara lo corresponden 3 puntos y a la cruz 4 puntos. En este 


caso. tenemos 
МӨ = 0,5-8 + 0,5-4 = 8,5; 
DO = 0,5 (32 + 42) — (8,5) = 0,25, 


Como vemos МӨ == Mx mientras que DO < Dx. Lo último era de 
esperar ya que los valores de @ puedon diferir de 8,5 en + 0,5 todo lo 
más mientras que la dispersión de los valores do x puede alcanzar 


hasta + 2,5. 


2,2. Variables aleatorias continuas. Supongamos que 
en el origen de coordenadas del plano se encuentra una can- 


1) Supongamos que además de la variable E estamos 
observando otra variable aleatoria 1. Si la distribución de la variable E 
no se altera al conocerse el valor que ha tomado la variable 4, es 
lógico decir que & no depende de 1). 
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tidad determinada de radio. Todo átomo de radio emite 
al desintegrarse wna partícula a. Su dirección puede ser 
caracterizada mediante el ángulo % (fig. 3). Tanto desde 


Fig з 


el punto de vista teórico como práctico, la emisión es po 
sible en cualquier dirección; por eso, esta variable alea- 
toria puede tomar cualquier valor comprendido entre 0 
y 2n. 

Diremos que una variable aleatoria Ё es continua si puede 
tomar cualquier valor comprendido en um intervalo (a, b). 

Toda variable aleatoria continua Ё queda definida si 
se da el inlervalo (а, b) que contiene los valores posibles 
de esta variable y la función p (x) que lleva el nombre de 
densidad de probabilidad de la variable aleatoria Ё (o den- 
sidad de distribución de E). 

El significado real de la función р (х) es el siguiente: 
sea (а”, b") un intervalo cualquiera contenido en (a, Б) (es 
decir, sea аса” y b' < b); entonces la probabilidad de 
que Ё tome un valor perteneciente al intervalo (a”, b') 
es igual a la integral 

br 
Pa Б) | рді. (44) 


“ 
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El área sombreada en la fig. 4 es igual al valor de la inte- 
gral (14). 

El conjunto de los valores de Ё puede formar un intervalo 
cualquiera. Puede darse incluso el caso de que а = — оо 


- 


о a е b' b 


хү 


a 
Fig. 4 


y también el de que b = оо. En cambio, la densidad р (=) 
debe cumplir dos condiciones análogas a las condiciones (1) 
y (2) para las variables discretas: 

a) la densidad р (т) debe ser positiva: 


р (ж) > 0; (15) 


b) 1а integral de la densidad р (x) correspondiente a todo 
el intervalo (а, b) debe ser igual a 1: 


b 
Voq dat. (16) 


Se denomina esperanza matemática de la variable aleato- 
ria continua el número 


b 
Mi- f xp (2) de, an 


El significado de esta característica өз el mismo que en 
el caso de una variable aleatoria discreta. Efectivamente, 
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puesto que " > 
M= { 2p(z) dz] | р (2) dz, 


es fácil ver que se trata del valor medio de Ẹ pues cualquier 
valor de Ё es un número т del intervalo (a, 5) que aparece 
en la integral con el peso p (x) 1). 

Todo lo expuesto en el punto 2.1, desde la fórmula (4) 
hasta la fórmula (13) inclusive, sigue siendo válido también 
para las variables aleatorias continuas; esto se refiere a la 
definición (7) de la varianza, a la fórmula (9) que se emplea 
para calcularla y a todas las propiedades de МЕ y ОЕ. La 
repiticién huelga. 

Mencionemos únicamente otra fórmula más que permite 
calcular la esperanza matemática de una función de E. 
Sea, como antes, p (z) la densidad de probabilidad de la 
variable aleatoria E. Tomemos una función continua cual- 
quiera f (г) y consideremos la variable aleatoria т =f (5). 
Se puede demostrar que 


b 
MfG = | рб) da. (18) 


Cabe subrayar que, hablando en términos generales, 
Mi (2) зе/ (ME). 

Si la variable aleatoría y está definida en el intervalo 
(0, 1) y su densidad es p (2) = 1 (fig. 5), se dice que ostá 
uniformemente distribuida en (0, 1). 

Efectivamente, cualquiera que sea el intervalo (а”, b') 
que se tome dentro de (0, 1), Ja probabilidad de que * tome 
un valor perteneciente al intervalo (a', b') es igual a 


” 
POLE 


1) También en este caso podemos señalar una fór- 
mula apáloga de la Mecânica: si p G) es la densidad lineal de la barra 
a 5 = < b, la abscisa del centro de gravedad de la barra se calcula 
mediante la fórmula 


z 


І 
ane 


xp (z) az] | p (z) dz. 
° 
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о sea, es igual a la longitud de este intervalo. En particular, 
si dividimos (0, 1) en un nümero arbitrario de intervalos 


Р(х) 


Fig. 5 


de igual longitud, la probabilidad de que y tome un valor 
perteneciente a cualquiera de estos intervalos será la misma. 
Es fácil calcular que 


i 1 
My= | zp(x)dz— | adam; 
0 0 


Dy= | 22р (2) dz — (My) Page 
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otm 


La variable aleatoria y será empleada frecuentemente 
өп lo sucesivo. 

2.3. Variables aleatorias normales. Se da el nombre de 
variable aleatoria normal (o gaussiana) a toda variable 
aleatoria { que está definida en todo el eje (—oo, оо) y que 
tiene la densidad 

po QUAE d 
r= aes (19) 


donde a y о > 0 son unos parámetros numéricos 1). 


1) а өз un número y no una variable aleatoria. Sin 
embargo, empleamos esta letra griega ya por tradición. 
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El parámetro а no influye en 1a forma de la curva p (2): 
su variación sólo conduce a un desplazamiento de la curva a 
lo largo del eje x. En cambio, al variar 0 sí se altera la forma 

de la curva. En efecto, es fácil ver que 

1 
máx = ————. 

р(а)--р(а) sy 
Si disminuimos 0, aumentará el máx р (z). Pero toda el 
Área que encierra la curva p (z) es igual a 1 segün (16). Por 


eso, la curva tendrá que estirarse hacia arriba en una vecin- 
dad de 2 = a y decrecer para todos los valores suficiente- 
mente grandes de т. En la fig. 6 hemos representado dos 
densidades normales correspondientes al caso a — 0 y 
c —1 y al caso a =0 y с = 0,5. (Otra densidad normal 
será construida en la fig. 21 de la página 54). 

Se puede demostrar que 


М5 =a y D$ =o. 


Las variables aleatorias normales aparecen con gran 
frecuencia al analizar cuestiones de la más diversa indole. 
La razón de esto será explicada a continuación. Por ejemplo, 
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el error de medición 8 es, como regla, una variable aleatoria 
normal. Si no existe error sistemático en las mediciones, se 
tieno а == Мӧ = 0. La magnitud c = V Dó, llamada 
error medio cuadrático, caracteriza el error del método que 
se emplea para la medición. 

La regla de «las tres sigmas». Es fácil calcular que 


24,36 
1 p (z) dz = 0,997 
2-30 
Cualesquiera que sean e y o en (19). De (14) resulta que 
Р {в — 80 < 5 < а + 30) = 0,997. (20) 


La probabilidad 0,997 es tan próxima a 1 que la ültima 
férmula se interpreta a veces asi: al efectuar una prueba es 
prácticamente imposible obtener un valor de G que difiera de 
ME en más de 3c. 

2.4. Teorema central del límite en la Teoría de las proba- 
bilidades. Laplace fue e] primero en enunciar este importante 
teorema. Varios matemáticos destacados, entre los cuales 
figuran P. L. Chébishev, A. A. Markov y А. M. Liapunov, 
estudiaron sus generalizaciones. Su demostración es bastante 
compleja. 

Consideremos № variables- aleatorias E,, En, ..., En 
independientes 6 idénticas de modo que las distribuciones de 
probabilidad de estas variables coinciden. Por consiguien- 
te, también coinciden sus esperanzas matemáticas y varian- 
таз. 
Pongamos 


МЕ, = ME, =... = МЕ = т, 
DE, = DE, =... = DEy = Б 
Designemos por Py la suma de todas estas variables: 
es =E tit... En 
De las fórmulas (6) y (13) se deduce que 
Mex = M (& + Es +... 5) = Nm, 
Dew = D (& + ba t... + En) = ND. 


Consideremos ahora la variable aleatoria normal Ey 
con estos mismos parámetros: а = Nm y o? = Nb? 
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El teorema central del limite afirma que cualquiera que 
sea el intervalo (a', b') se tiene para valores grandes de N 
M 
Р (o! «os <b} ле p, (2) de. 


El significado real de esto teorema өз obvio: la suma py 
de una gran cantidad de variables aleatorias idénticas es 
aproximadamente normal (рр, (2) 2 Pry (2). 

En realidad, este teorema es válido en situaciones mucho 
más generales: los sumandos Er, Б... Ey pueden no ser 
idénticos e independientes; lo que importa es sólo que cada 
sumando por separado no influya considerablemente en 1а 
suma. 

Es este teorema el que explica porque las variables alea- 
torias normales aparecen con tanta frecuencia en la práctica. 
En efecto, siempre que tropezamos con el efecto sumario de 
una gran cantidad de factores aleatorios insignificantes, la 
variable aleatoria resultante es normal, 

Por ejemplo, la desviación del proyectil respecto al obje- 
tivo es casi siempre una variable aleatoria normal por cuanto 
depende de las condiciones meteorológicas en las diferentes 
partes de la trayectoria y de otros muchos factores. 

2.5. Esquema general de aplicación del método de Monte- 
carlo. Supougamos que debemos calcular una magnitud m 
que desconocemos. Tratemos de idear una variable aleatoria 5 
tal que ME = т. Sea además DE = b°. 

Consideremos № variables aleatorias independientes É,, 
$4. En con la misma distribución que tiene Ё. Si 
N es Aufl ensem ante grande, la distribución de la suma 
Pw = E+ E+ . + Ey será, según el teorema del 
puits 2.4, una distribución aproximadamente normal con 

os parámetros а = Nm y о? = Nb*. De (20) se deduce que 


Р (Nm — 8bV N « py < Nm + 8bY М} ~ 0,997. 
Si dividimos entre N la desigualdad qne figura ontre las 
llaves, obtendremos una desigualdad equivalente; la proba- 
bilidad continuará siendo la misma: 
Ed E 
P (77 <= MUS) яе 0,997. 


3-1216 
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Podemos representar este Ultimo resultado en la forma si- 
guiente 


N 
P (lr тр) > 0,997. (21) 
jmi 


Ésta es una relación de suma importancia para el método 
de Montecarlo. Ofrece un método para calcular m y a le vez 
permite estimar el error. 

En efecto, determinemos N valores de la variable aleato- 
ria E ®). Según la fórmula (21), la media aritmética de estos 
valores será igual aproximadamente a m. Con una elevada 
probabilidad podemos afirmar que el error de esta aproxima- 
ción no pasa de 3b/V/ М. Es obvio que el error tiende a cero 
cuando N crece. 


#3. 
OBTENCION 

DE VARIABLES ALEATORIAS 
EN LAS MCE 


El propio planteamiento de la cuestión («obtención de 
variables aleatorias en las MC E») provoca cierta perplejidad: 
ide qué casualidad se puede hablar si todo cuanto realizà la 
máquina debe ser programado de antemano? 

En esta cuestión hay efectivamente ciertas sutilezas; 
pero conciernen más bien a la filosofía y, por eso, las deja- 
remos a uri lado. 

Sin embargo, cabe sefialar que las variables aleatorias 
consideradas en el $ 2 representan conceptos matemáticos 
ideales. Sólo la práctica dilucida si pueden ser empleados 
para describir determinados fenómenos de la naturaleza. 
Esta descripción es siempre aproximada. Es más, una va- 
tiable aleatoria que describe satisfactoriamente una magnitud 
física en unos fenómenos, puede resultar uma característica 


1) Resulta lo mismo determinar un valor de cada 
una de las variables £j, £j, . - ., Ex о hallar № valores de la variable E 
ya que todas estas variables aleatorias son idénticas (tienen la misma 
distribución). 
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desafortunada de esta misma magnitud al tratarse de otros 
fenómenos. 

De la misma forma, un camino que se asomeja a una 
recta (a una recta matemática ideal «sin anchura») en el 
mapa del país, se convierte en una franja sinuosa en el 
plano a gran escala de una localidad... 

Suelen distinguirse tres procedimientos que se aplican 
para obtener variables aleatorias: tablas de nümeros alea- 
torios, generadores de números aleatorios y el método do 
números seudoalentorios. 

3.1. Tablas de números aleatorios. Realicemos el өхрегі- 
mento siguiente. Tomemos diez tarjetas iguales escribiendo 
en ellas las cifras 0, 1, 2, . . ., 9, metamos todas las tarjetas 
en un sombrero y, después de mezclarlas bien, comencomos a 
extraerlas, con la particularidad de que cada vez sacaremos 
una tarjeta solamento, la volveremos a meter en el sombrero 
y mezclaremos de nuevo todas las tarjotas. Las cifras así 
obtenidas pueden ser recogidas en una tabla semejante a la 
tabla A de la página 78 al final del libro (por razones de 
comodidad las cifras en la tabla A vienen agrupadas de cinco 
en cinco). 

Esta tabla lleva el nombro de tabla de números aleatorios 
aunque sería más justo denomiuarla tabla de cifras aleato- 
rias. Dicha tabla puede ser introducida en la memoria de la 
MCE, Cada vez que, durante los cálculos, necesitemos un 
valor de la variable aleatoria con la distribución 


0 4 grmi 
hu 04 04 ad (22) 


podemos tomar una tras otra las cifras de esta tabla. 

La mayor tabla de nümeros aleatorios publicada hasta el 
momento contiene 1 000 000 de cifras!). Por supuesto, fue 
obtenida empleando, en lugar del sombrero, una técnica más 
moderna: una ruleta electrónica. Un esquema elemental de 
esta ruleta aparece en la fig. 7 (el disco giratorio es para- 
do en seco y se toma la cifra que señala la flecha). 

Es preciso subrayar que no es tan fácil como parece obte- 
ner una tabla buena de números aleatorios, Todo aparato 


^ 1) RAND Corporation, A million random digits 
with 100 000 normal deviates, The Free Press, 1955. 


” 
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físico real genera variables aleatorias cuyas distribuciones 
difieren ligeramente de la distribución ideal (22). Además, 
durante el experimento pueden surgir errores (por ejemplo, 
una de las tarjetas puede adherirse. por cierto tiempo al 
forro). Por eso, las tablas obtenidas suelen comprobarse 
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Fig, 7 


minuciosamente a base de tests estadísticos especiales para 
ver si no hay contradicción entre unas u otras propiedades 
de grupos de números y la hipótesis de que estos números 
representan valores de la variable aleatoria (22). 

Veamos uno de los tests más sencillos. Consideremos una 
tabla qua contiene М cifras. Supongamos que el número de 
ceros en esta tabla es v, el número de unos es v,, el número 
de doses es va, ete. Calculemos la suma 


9 
2 (9, —0,4 NY. 
(о 

La Teoría de las probabilidades permite predecir los 
límites entre los cuales puede estar comprendida esta suma 
que no puede ser excesivamente grande (por cuanto la espe- 
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тапта matemática de cada v, es igual a 0,1 ^) ni excesiva- 
mente pequefia (ya que esto equivaldria a una «regularidad 
excesiva» en la distribución de los valores). 

Las tablas de nümeros aleatorios se emplean sólo cuando 
les cáleulos correspondientes al método de Montecarlo se 
realizan a mano. Esto se debe а que todas las MCE poseen una 
memoria interna relativamente reducida en la que no se 
puede introducir una tabla grande. Por otra parte, tampoco 
conviene introducir la tabla en la memoría externa y recu- 
rrir constantemente a ella ya que esto afectaría considerable- 
mente la velocidad del cálculo. 

Sin embargo, no está excluido que, con el tiempo, la 
capacidad de memoría de las MCE aumente sensiblemente 
y entonces las tablas de nümeros aleatorios resulten de gran 
utilidad. 

3.2. Generadores de números aleatorios. La ruleta mencio- 
nada en el punto 3,1 podría, al parecer, combinarse con la 
computadora y emplearse para generar, según sea necesario, 
números aleatorios. Pero todo dispositivo mecánico resulta 
demasiado lento comparado con la MCE. Por eso, para 
generar variables aleatorias suelen emplearse los ruidos 
en las lámparas electrónicas: si durante un intervalo fijo At 
de tiempo el nivel del ruido sobrepasa el umbral escogido 
un número par de veces, se inscribe el cero; si esto ocurre 
un número impar de veces, se inscribe el uno!), 

A primera vista, éste es un procedimiento muy cómodo. 
Imaginemos m generadores de este tipo que funcionan sin ce- 
sar paralelamente y envían ceros y unos aleatorios a todos 
los órdenes binarios de una célula especial. Ceda ciclo co- 
rresponderá а un número de m órdenes. Durante el cálculo 
podemos dirigirnos en cnalquier momento a esta célula y 
tomar de ella un valor de la variable aleatoria y uniforme- 
mente distribuida en el intervalo (0, 1), Este valor, claro 
está, será aproximado y quedará representsdo por una frac- 
ción binaria de m órdenes 0, «ау - · . Gem, donde cada 
Ош imita la variable aleatoria con la distribución 


| о 1 
Ж 


3) Existen procedi 


entos más perfectos, 
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Sin embargo, también este método tiene sus defectos. 
En primer lugar, es dificil comprobar Ja «calidad» de los 
nümeros generados. Se hace necesario un control periódico 
ya que debido a distintos desarreglos puede surgir la así 
llamada «deriva en la distribución» (cuando los ceros y los 
unos en uno de los órdenes dejan de aparecer con 1a misma 
frecuencia). En segundo lugar, todos los cálculos en las 
MCE suelen realizarse dos veces para evitar intermitencias 
casuales. Pero, sin recurrir durante los cálculos a la memo- 
ria, es imposible reproducir los mismos nümeros aleatorios 
y si recurrimos a la memoria, volvemos al caso de las tablas. 

Generadores de este tipo resultarán itiles, sin duda, 
cuando comenzarán a construirso MCE especializadas, adap- 
tadas a la solución de problemas por el método de Monte- 
carlo. En las MCE universales, que raramente se aplican 
para los cálculos a base de números aleatorios, no es econó- 
mico mantener y explotar un dispositivo especial. Resulta 
más conveniente emplear los así llamados números seudoalea- 
torios. 

3.3. Nümeros seudoaleatorios. Como quiera que la «cali- 
dad» de los nümeros aleatorios empleados se controla median- 
te tests especiales, podemos no interesarnos por el procedi- 
miento aplicado para obtenerlos exigiendo solamente que 
verifiquen el sistema de tests adoptado. Incluso podemos 
intentar calcularlos a partir de una fórmula dada. Por su- 
puesto, debe ser una fórmula muy ingeniosa. 

Los nümeros que se obtienen a partir de una fórmula y que 
imitan los valóres de la variable alestoria y se denominan 
nümeros seudoaleatorios. La palabra «imitan» significa que 
estos números verifican una serie de tests como si fuesen 
valores de esta variable aleatoria. 

El primer algoritmo destinado a la construcción de nüme- 
ros seudoaleatorios fue propuesto por J. von Neumann y se 
conoee como el método de centros de cuadrados. Para explicarlo 
recurriremos a un ejemplo. 

Consideremos el número y, = 0,9876 formado por cuatro 
cifras. Al elevarlo al cuadrado obtendremos el número 
үз = 0,97535376 de ocho cifras. Tomemos las cuatro cifras 
que aparecen en el centro de este número y pongamos y, = 
= 0,5353. 

Elevando ahora y, al cuadrado (y? = 0,28654609), escoja- 
mos de nuevo las cuatro cifras del centro y pongamos y, = 
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х= 0,6546. Procediendo do la misma forma, tendremos y? = 
= 0742850116, y, = 0,8501; y = 0,72207001, y, = 0,2670; 
yi = 0,07128900, ys = 0,1289, etc.!). 

Pero este método no puede considerarse satisfactorio pues 
da más números pequeños de lo que se necesita, Por eso, 
fueron desarrollados por distintos científicos otros algoritmos. 
Algunos se basan en las particularidados de las MCE con» 
cretas. A título de ejemplo, consideremos uno de estos algo- 
ritmos empleado en la MCE Strelá. 


Ejemplo, 2). Strelá es una M CE de tres direcciones con coma flolan- 


te, La célula en la que se inscribe el número z consta de 43 órdenes bina- 
rios (fig. 8). La computadora opera con nümeros binarios en forma 


normalizada: z = + 4.2*?, donde p es el orden y 4 la mantisa del 


o T2 T4 TS] Үз: 32[33[5а]55[з36[37]38]59[45]а1]52 


Montiso } Orden 


Signo de montisa Signo del orden 


Fig. 8 


número. En el orden j-ósimo de la célula puede figurar un cero o un 
uno; representemos por в; esta variable. Entonces 


121 Li 535 
«етене Es Y 
D = egr?" -H 635224... A £520. 


La condición de normalización 0,5 < q < 1 significa que e; es siem- 
pre igual a 1. El signo « + > se representa por el cero y el signo <-> 
por el uno. 
_ El número үү se obtiene del número y, realizando tres opera- 
ciones: 

1) yy se multiplica por una constanto grande que, como regla, 
so toma igual a 10%; 


1) Podemos resumir este algoritmo en la fórmula 
Yapi = Fn). donde F representa el conjuuto de operaciones que 
debemos realizar con el número y, para obtener урау. El número yg eS 
un número dado. 

2) Véase И. М. Соболь, Псевдослучайные числа для машины 
Cmpesa, Теория веролтности и ве примепешия, 8, Ni 2 (1958), 205— 
214 (I. M. Sóbol, Números soudoaleatorios para lu computadora Strelá; 
en la revista «Teoría de las probabilidades y sus aplicaciones», 3, 
n° 2 (1958), pág. de 205 a 241). 
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2) la imagen del producto 1017 үу se desplaza en siete órdenes hacia 
la izquierda (de modo que los sieto primeros órdenes del producto desa- 
parecen miontras que en los órdenes del 36 al 42 aparecen ceros); 

3) se toma el valor absoluto del número obtenido (después de nor- 
marlo) que será precisamente Үлы. 


À arrancamos de yo = 1, este proceso permite obtener más de 
80 000 números diferentes yy; después en la sucesión surge un período 
y los números comienzan a repetirse, Diferentes controles de los 
50 000 primeros números han dado resultados plenamento satisfacto- 
rios, Estos números вө hau empleado más de una vez para la solución 
de diversos problemas, 


Las ventajas del método de números seudoaleatorios 
son evidentes. En primer lugar, para obtener un número 
basta realizar unas cuantas operaciones sencillas y gracias 
a esto la velocidad de generación de números aleatorios es 
del mismo orden que la velocidad de funcionamiento de la 
МСЕ. En segundo lugar, el programa ocupa sólo unas cuan- 
tas células de la memoria, En tercer lugar, cualquiera de 
los números y, puede ser reproducido fácilmente. Por últi- 
mo, basta comprobar una vez solamente la «calidad» de esta 
sucesión para poder aplicarla después reiteradamente y con 
seguridad en la solución de problemas semejantes. 

El método tiene sólo un defecto: la «reserva» de los nú- 
meros seudoaleatorios es limitada. Pero existen procedi- 
mientos que permiten obtener una cantidad mucho mayor 
de números. En particular, se pueden variar los números 
iniciales yo. Actualmente los cálculos en los que se emplea 
el método de Montecarlo se realizan, en su aplastante mayo- 
ría, aprovechando los números seudoaleatorios. (Véanse in- 
formaciones complementarias en le $ 10). 


$4. 
"TRANSFORMACIONES i 
DE VARIABLES ALEATORIAS 


En la solución de diferentes problemas tropezamos con la 
nocesidad do simular distintas variables aleatorias. En las 
etapas iniciales del desarrollo del método de Montecarlo los 
investigadores intentaban obtener cada variable aleatoria 
construyendo una ruleta apropiada. Por ejemplo, para hallar 
los valores de la variable aleatoria con la distribución 


siguiente 


“4 от, їз д , 
los 0,25 0,125 oss) (22) 


se puede emplear la ruleta representada en la fig. 9 que fun- 
ciona igual que la ruleta de la fig. 7, pero no está dividida 
uniformemente sino en partes proporcionales a las py. 


Fig. 9 


Pero esto no hace falta: resulta que los valores de cual- 
quier variable aleatoria se pueden obtener efectuando trans- 
formaciones de una sola variable aleatoria (es decir, de una 
variable estandard»). Suele emplearse con este fin la variable 
aleatoria y uniformemente distribuida en el intervalo (0, 1). 
Conocemos ya cómo se obtienen los valores de y. 

Denominaremos sorteo de la variable aleatoria E ol pro- 
ceso de determinación de un valor de esta variable aleatoria Ё 
mediante la transformación de uno o varios valores de y. 

4.1. Sorieo de una variable aleatoria discreta. Suponga- 
mos que es preciso obtener valores de la variable aleatoria E 


con la distribucién 


& n ы Құда 
Тір р... ра 


Consideremos el intervalo 0 < y < 1 dividiéndolo en n 
intervalos de longitudes p;, Pa, .. ., рл. Es evidente que 
las coordenadas de los puntos de isión serán y = p; 
Y =P +P» У = р раі Pye. YS Pt Pet 
+... + pa-; Indiquemos con las cifras 1, 2,,.., n los 


> 
9 PL рер, Py Py Ps EA Vy 
Fig, 10 


intervalos obtenidos (fig. 10). Con esto concluyen los prepa- 
rativos del sorteo de E. Cada vez que tengamos que «reali- 
zar un experimento» y sortear el valor de E, tomaremos un 
valor de y y construiremos el punto y = y. Si este punto 
aparece en el intervalo correspondiente al número i, acep- 
taremos que E = z; (en este experimento). 

Es muy fácil argumentar la legitimidad de este procedi- 
miento. En efecto, puesto que la variable aleatoria y está 
uniformemente distribuida en (0, 4), la probabilidad de que 
y pertenezca a uno de los intervalos es igual a Ja longitud 
del mismo. Por consiguiente, 


P(0cv«cp)—p. 


Segün nuestro método, tomamos ё = 2; cuando 
Pi тра +... Бра <р + ра +... рг 


y la probabilidad de ello es igual a pi. 

Por supuesto, al emplear la MCE podemos pasarnos sin 
la fig. 10. Supongamos que tanto los námeros Zy, zs, ..., Ln 
como las probabilidades py, p, -+ Poy р; + pa + pss +...» 1 
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han sido introducidos uno tras otro en la memoria. En la 
fig. 11 de la pág. siguiente representamos e! esquema 
sinóptico del subprograma de sorteo de E. 


Ejemplo. Sortéense 10 valores de la variable aleatoria con la dis- 
tribución 


з 4 
8 (5% 0,42) * 

Escojamos los valores de y tomando pares de cifras de la tabla A 
(pág. 78) multiplicados por 0,01 !). Es decir, y = 0,86; 0,51; 0,59; 
0,07; 0,05; 0,66; 0,45; 0,56; 0,04; 0,34. 

Según nuestro esquema, los valores de y menores que 0,58 corres- 
ponden a 0 = 3 mientrus que los valores y > 0,58 corresponden a 0 = 


= 4; de modo gue obtenemos los valores siguientes: 0 = 4) 33 4; 3; 4; 
415131418. 

Cahe subrayar que el orden de numeración de los valores 
ац, Eg, - +, 2, en la distribución de § puede ser arbitrario 
pero debe fijarse antes de comenzar el sorteo. 

4.2. Sorteo de una variable aleatoria continua. Suponga- 
mos ahora que debemos obtener valores de una variable 
aleatoria & distribuida con la densidad p (z) en el interva- 
lo (a, 5). 

Demostremos que los valores de E so pueden determinar 
de la ecuación 


| р(х) dx ==, (23) 


a 


o sea, que escogido ya el valor do y es preciso revolver la 
ecuación (23) para hallar el valor de E. 
Para demostrar esta afirmación consideremos la función 


y= joe dz. 


1) Puesto que on este pals los valores de p, tienen 
dos cifras decimales, basta tomar los valores de y con dos cifras 
decimales. Aceptado este grado de aproximación, so nos puede dar el 
caso de y = 0,58 que debe ser agregado al caso de y > 0,58 (ya que 
el valor y= 0,00 es posible mientras que el valor y = 1,00 es 
imposible). Si empleamos para y números con varias cifras decimales, ol 
сазо ү —p, se hace poca probable y puede ser agregado a cualquiera 
de las desigualdades. 
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Determinamos y 


Tomamos 
ізт 


Aumentamos en 1 las direc- 
ciones de las células de las 
que se toman py y 21 


Reconstruimos las 
direcciones 


de руу х4 


Fig. 14 
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De las propiedades generales (15) y (16) do la densidad resulta 
quo 


y (a) = 0, y (b) —1 
y que la derivada 
у (2) = p (2) >0. 


Por lo tanto, la función y (x) crece monótonamente del O 
al 1 (fig. 12). Debido a esto, toda recta y = y, donde 0 << 
< ү << 1, corta е) gráfico de y = y (1) on un punto único 


1 


1 


Fig. 12 


cuya abscisa se toma precisamente igual a E. Es decir, la 
ecuación (23) tiene siempre solución y esta solución es ünica. 

Tomemos ahora un intervalo arbitrario (a', b^) conte- 
nido en (a, b). A los puntos 


а <2< 6 


de este intervalo les corresponden las ordenadas de la curva 
y = у (z) que satisfacen la desigualdad 


y (a^) < y « y (b^). 


Por eso sí E pertenece al intervalo а” < z < b', resulta 
que y pertenece al intervalo y (a^) « у < y (b') y viceversa 


(fig. 13). Luego 
P< Sch} =P (y (a) <y< y (6). 


Puesto que y está unifermemente distribuida en (0, 1), 
lenemos 


v 
P (Gn) ca 00 a) 7) — | p (a) dz. 
Es decir, 
ы 


Ра <k<ly= | (аш; 


pero esto significa precisamente que la variable aleatoria 


9 о ав ы” 


Fig. 13 


E, raíz de la ecuación (23), tiene la densidad de probabilidad 
р (3). 
Ejemplo. Se dice que la variable aleatoria y está uniformemente dis- 


инал en el intervalo (a, b) si su densidad es constante en este inter- 
valo: 


ғ()=р рага acc, 
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Para sortear los valores de т formamos la ecuación (23): 


do donde obtenemos una expresión explícita para m: 


n=0+ (6 a). (24) 


Otros ejemplos de aplieación de la fórmula (23) pueden 
verse en los puntos 5.2 y 8.3. 

4.3. Método de Neumann para el sorteo de una variable 
aleatoria continua. Puede resultar que sea muy difícil resol- 
ver la ecuación (23) respecto a E; por ejemplo, si la integral 
de p (z) no se expresa en funciones elementales o si la den- 
sidad p (z) viene dada gráficamente. 

Supongamos que la variable aleatoria E está definida 
en un intervalo finito (a, b) y que su densidad está acotada 
(fig. 14): 


p(nsM,. 


Los valores de 2 se pueden sortear de la siguiente manera 

1) se escogen dos valores y” y y” de la variable aleatoria y 
y se considera el punto aleatorio Г (n’, 1”) con las coor- 
denadas 


m =a t+ y (b—a)ym'—y"My 


2) si el punto Г aparece por debajo de la curva y = 
= p (т), tomamos Ё = n’; si el punto Г aparece por encima 
de la curva y = р (г), retiramos de la consideración el 
par (y’, y”) y escogemos un par nuevo de valores (y', y^^. 

La argumentación de este método viene en el punto 9.1. 

4,4, Sobre el sorteo de variables normales. Existe una 
gran variedad de procedimientos que se aplican para sortear 
distintas variables aleatorias. No nos detendremos en ello. 
Suelen emplearse sólo si los mótodos de los puntos 4.2 y 4.3 
resultan poco efectivos. 
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En particular, esto ocurre en el caso de la variable alea- 
toria normal £ ya que la ecuación 


5 
va) 
VE 
по admite solución explícita y el intervalo de los valores posi- 
bles de б es infinito. 


E 
Стй-ү 


YEP (x) 


«Y 


Fig. 14 


La tabla B que viene al final del libro (en la pág. 78) 
contiene valores (ya sorteados) de la variable aleatoria nor- 
mal ¢ con la esperanza matemática МЁ = 0 y la varianza 
05 = 1. Es fácil demostrar!) que la variable aleatoria 


б =a + оё (25) 
también será normal; además, de (10) у (11) se deduce que 
Mo’ =a y DY = а? 


1) La demostración viene en ol punto 9.2. 
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Por consiguiente, la fórmula (25) permite, empleando la 
tabla B, sortear valores de cualesquiera varinbles normales. 

4.5. Otra vez sobre el ejemplo del punto 1.2. Ahora pode- 
mos explicar cómo se han escogido los puntos aleatorios en 
las figuras 1 y 2. En la fig. 1 hemos construido los puntos 
con las coordenadas 

к= egy" 

Los valores y” y Т” se han calculado a base de los grupos 
de cinco cifras de la tabla А: х, = 0,86515, y, = 0,90795; 
Ty = 0,66155, y, = 0,668434, otc. 

Se puede demostrar!) que, debido a que las abscisas y las 
ordenadas de nuestros puntos son independientes, la probabi- 
lidad de que uno de estos puntos aparezca en una región 
cualquiera interior del cuadrado es igual al área de esta 
región. En otras palabras, los puntos están uniformemente 
distribuidos en el cuadrado. 

En la fig. 2 hemos construido los puntos con las coordens- 
das 

z =0,5 + 0,26" e y =0,5 + 0,207" 
con la particularidad de que los valores {' y £'* se toman 
de la tabla B uno tras otro: 
z, = 0,5 + 0,2-0,2005, y, = 0,5 + 0,2 -1,1922; 
x, = 0,5 + 0,2-(—0,0077), ... 
Un punto no se ha considerado por quedar fuera del cuadra- 
do. 
De la fórmula (25) se deduce que las abscisas y las ordena- 


das de estos puntos representan variables aleatorias normales 
con la esperanza matemática а = 0,5 y la varianza a* = 0,04 


* La demostración viene en el punto 9,3. 


CAPITULO 2 
EJEMPLOS 
DE APLICACION 
DEL MÉTODO 
DE MONTECARLO 


$8. 
ANALISIS 
DE UN SISTEMA DE SERVICIOS 


5.1. Descripción del problema. Consideremos uno de 
los sistemas más sencillos de servicios. Consta de m líneas 
(canales o centros de servicio), cada una de las cuales se 
puede emplear para «atender a los consumidores». El sistema 
recibe pedidos, con la particularidad de que los momentos 
de recepción son aleatorios. Todo pedido va a la línea núme- 
ro 1. Sien el momento de recepción del l-ésimo pedido (de- 
signemos este momento por 7,) la línea está libre, comienza 
a cumplir el pedido durante un tiempo Zoe (fo, es el tiempo de 
ocupación de la línea). Si en el momento F, la línea número 1 
está ocupada, el pedido se transmite instantáneamente a la 
línea número 2, etc. Por último, si en 61 momento Т, resul- 
tan ocupadas todas las líneas, el sistema rechaza el pedido. 

Es preciso determinar el promedio de pedidos que atenderá 
y rechazará el sistema durante un tiempo 7. 

Está claro que semejantes problemas pueden surgir al 
analizar el funcionamiento de cualquier empresa y no sólo 
de las empresas de servicios. En algunos casos muy espe- 
ciales se logra la solución analítica. Pero, en casos complejos 
(de los que trataremos más tarde), el método de Montecarlo 
constituye el único método de cálculo. 

5.2. Flujo elemental de pedidos. El primer problema que 
se plantea al analizar sistemas de este tipo es: ¿qué repre- 
senta el flujo de pedidos? Este problema se resuelve de mane- 
ra experimental a base de un estudio suficientemente pro- 
longado de los pedidos. El análisis de los flujos de pedidos, 
correspondientes a distintas condiciones, ha permitido desta- 
car algunos casos que suelen darse con suficiente frecuencia. 

Un flujo de pedidos se denomina flujo elemental (o flujo 
de Poisson) si el intervalo de tiompo т que transcurre entre 
dos pedidos sucesivos representa una variable aleatoria 
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distribuida en el intervalo (0, oc) con la densidad 
p (x) =ae-%, (26) 
La ley (26) lleva también el nombre do distribución ezponen- 


cial (véase la fig. 15, donde hemos representado las densi- 
dades (26) para а — 1 y a = 2). 


о=2 


Fig. 15 


Es fácil caleular la esperanza matemática de v: 


Mt= f гр (2) dr= { тае-°® dz. 
D 5 


Integrando por partes (и = x y dv = ae-**dz), obtenemos 


= 
=ar 
Mt— [2061 | ent dz [ -———h-i 
o 
El parámetro a se denomina densidad del flujo de pedidos. 
La tórmula correspondiente al sorteo de т se obtiene 
fácilmente de la ecuación (23) que tiene en nuestro caso la 
forma 
+ 
f ae~ dx =y. 
% 
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Calculando la integral del primer miembro, obtenemos la 
relación 
d—e-—, 
de donde resulta 
_ 1 
t= —.h0-». 


Como quiera que la variable 1 — ү tiene la misma distri- 
bución que y, podemos emplear en lugar de la última fórmula 
esta otra 


t= —tiny. (27) 


5.3. Esquema de cáleulo. Consideremos, pues, el fun- 
cionamiento del sistema descrito en el punto 5.1 en el caso 
de un flujo elemental de pedidos. 

Pongamos cada linea en correspondencia con una célula 
de la memoria interior de la MCE inscribiendo en la célula 
el momento en que queda libre la linea. Sea %; el momento 
en que queda libre la i-ésima línea. Aceptaremos que los 
cálculos se inician en el momento 7, = 0 en el que se recibe 
el primer pedido. Podemos aceptar que en este momento 
todos los / son iguales a Т, ya que todas las líneas están 
libres. El tiempo final del cálculo será Ty, = 7, + T. 

El primer podido va a la línea número 1. Por consiguien- 
te, esta línea permanecerá ocupada durante el tiempo toc- 
Por eso, debemos tomar en lugar de 2 el valor nuevo (й) у = 
= T; + toc, agregar uno al contador de los pedidos cum- 
plidos y pasar a considerar el segundo pedido. 

Supongamos que hemos considerado уа k pedidos. De- 
bemos entonces sortear el momento de recepción del (k + 1)- 
ésimo pedido. Con este fin escogemos el valor siguiente de y 
y calculamos, a partir de la fórmula (27), el valor siguien- 
tet = т; después determinamos el momento de recepción 


Tya = Tyr. 


¿Estará libre para este momento la primera línea? Para 
responder, debemos comprobar la condición 


Ты. (28) 
Si esta condición se cumple, ello significa que en el momento 
Тұ la línea está ya libre y puede cumplir el pedido. Debe- 
remos sustituir й por Тун) + too, agregar uno al contador 
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de pedidos cumplidos y pasar a considerar el pedido si- 
guiente. 

Si 1а condición (28) no se cumple, ello significa que la 
primera linea está ocupada en el momento Ty. Comproba- 
mos entonces si está o no libre la segunda línea: 


fj «Ta. (29) 


Бі se cumple la condición (29), sustituimos /, por 7,44 + 
+ toc, agregamos uno al contador de pedidos cumplidos 
y pasamos a considerar el pedido siguiente. 

Si la condición (29) no se cumple, comprobamos la condi- 
cién 


tss Tu 
Puede ocurrir que para todo ¿ del 1 al n se tenga 


u> Tau 


o sea, puede suceder que todas las líneas estén ocupadas en 
el momento 7,4,. En este caso, agregamos uno al contador 
de pedidos rechazados y pasamos a considerar el pedido 
siguiente. 

Cada vez que caleulemos 7,,,, deberemos comprobar 
además la condición de finalización del experimento 


Tia > Trin 


Si esta condición se cumple, finalizamos el experimento. 
En los contadoras de pedidos cumplidos y de pedidos recha- 
zados tendremos entonces los números Beum Y Шгеев- 
Este experimento se repite М veces (empleando distintos 
valores de y) y se toma la media de los resultados encontrados 


N 
1 
Mecum >2-- У нањо» 
жа 


N 
1 
Mb rec FS -- Dy Breon d» 
jet 
donde cum h Y Heecn «p. Tepresentan los valoras de Leu 
y de п, en el j-ésimo experimento. 
En la fig. 16 aparece el esquema sinóptico del programa 
correspondiente a este cálculo (si es necesario, podemos obte- 
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Entrada de dato | 


b 


| іі 


i 


| Т0 type tg. Ht m0; kd | 
no 
Fin del experimento |e-—— ———| ThT pn? = ——. 


1 


ау = jj d | 


| 


E ns < № 
sí 
no 
(ny m Т ос | | db? 
{ по 
y 


[ 
"| demit 
y 


Agregamos 1 al contador 
de pedidos rechazados 


Agregamos 1 al contador de 
pedidos cumplidos 


; H 4 
n= 1 Iny 


Cálculo de resultados, salida 
de resultados, [in del 
problema Tra m Th. 
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ner en la sección «Fin del experimento» los valores de pe um (p 
y de uec i; que resultan en cada uno de los experimentos). 

5.4. Problemas más complejos. Es fácil persuadirse de 
que este mismo método se puede aplicar en sistemas mucho 
más complejos. 

Por ejemplo, £4, puede ser una variable aleatoria, y no 
una magnitud constante, y tomar diferentes valores para dis- 
tintas líneas (lo que significa que las líneas se distinguen 
por sus equipos o por la calificación del personal). El es- 
quema de cálculo, en esencia, continuará siendo el mismo, 
pero habrá que sortear cada vez los valores de to, correspon- 
diendo a cada línea una fórmula de sorteo propia. 

También se pueden considerar los llamados sistemas 
con espera que, en lugar de rechazar inmediatamente el 
pedido que no pueden cumplir, lo conservan durante cierto 
tiempo tp (tiempo de permanencia del pedido en el sistema); 
si en el transcurso de este tiempo se libera una de las líneas, 
dicha línea atiende el pedido. 

Además se pueden considerar sistemas оп los que el 
pedido en turno es atendido por la línea que se libera antes 
que otras. Por otra parte, es factible lomar en consideración 
las intermitencias aloatorias de las líneas, el tiempo aleato- 
rio de reparación de las mismas, la variación según el tiempo 
de la densidad del fluja de pedidos y otros muchos facto- 
res. 

Por supuesto, el análisis de estos sistemas exige deter- 
minados esfuerzos. Para obtener resultados de valor práctico 
debe ser ideado un modelo correcto. Esto requiere un estudio 
suficientemente minucioso de los flujos reales de pedidos, el 
cronometraje del funcionamiento de las secciones del sistema, 
otc. 

Hablando en términos generales, es preciso conocer las 
leyes probabilísticas de tuncionamiento de las secciones del 
sistema. En este caso, el método de Montecarlo permite ha- 
llar las leyes probabilísticas del funcionamiento de todo el 
sistema por complejo que sea. 

Estos métodos de análisis resulten muy ütiles cuando se 
trata de proyectar una empresa: en lugar de realizar un expe- 
rimento real costoso (y a veces imposible de realizar), pode- 
mosexporimentar en la MCE simulando distintas varian- 
tes de organización del trabajo o del empleo de los 
equipos. 
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8&0. 
ANALISIS 
DE LA CALIDAD 
Y DE LA SEGURIDAD 
DE PIEZAS 


6.1. Esquema elemental de análisis de la calidad. Consi- 
deremos una pieza 8 que consta de cierta cantidad (posible- 
mente grande) de elementos. Por ejemplo, si S es un aparato 
eléctrico, sus elementos pueden ser las resistencias (Боз), 
las capacidades (Cy), las lámparas, etc. Supongamos que 


22 К25%2У/ 


Fig. 47 


la calidad se determina por un parámetro 17 que puede ser 
calculado si se conocen los parámetros de todos los elementos 


U =] (Ran Fas 3 Са» Cas s de (0) 


Por ejemplo, si Ves la tensión en una sección de un circui- 
to eléctrico, podemos considerar, basándonos en la ley de 
Ohm, las ecuaciones del circuito y determinar U de estas 
ecuaciones. 

Sin embargo, los valores reales de los parámetros de los 
elementos difieren de los valores nominales. Por ejemplo, la 
resistencia representada en la fig. 17 puede tomar un valor 
cualquiera comprendido entre 20,9 y 23,1 kQ. 

Surge la cuestión siguiente; ¿cómo influyen en el valor 
de U las desviaciones que se observan entre los valores 
reales y nominales de los parámetros de todos los ele- 
mentos? 

Los límites de variación de U podrían estimarse escogien- 
do para todo elemento los valores «peores» de los parámetros. 
Pero no siempre ni mucho menos se conoce qué colección 
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de parámeiros será la «реог». Además, habiendo muchos 
elementos, esta estimación resultará exagerada por cuanto 
өз poco probable que todos los parámetros a la vez resulten 
«malos», 

Por eso, es más razonable tratar de estimar la esperanza 
metemática MU y la varianza DU considerando que tanto 
los parámetros de los elementos como el propio parámetro Y 
sou variables aleatorias. La magnitud M U representará el 
valor medio de U en toda la partida de piezas y la magnitud 
DU indicará qué desviaciones entre U y MU pueden pre- 
sentarse realmente. 

Recordemos que (según hemos señalado en el punto 2.2) 


MU 7 (MRa» MRs «+ 3 МСа» МС... od 


Siendo f una función más o menos compleja, se hace impo- 
sible el cálculo analítico de la distribución de U, A veces 
se logra determinarla experimentalmente analizando una 
gran partida de piezas elaboradas. Pero esto no siempre es 
posible y jamás puede realizarse en la fase de diseño de una 
pieza. 

Intentemos aplicar el método de Montecarlo. Para ello es 
necesario; a) conocer las características probabilísticas de 
todos los elementos y b) conocer la función f (más exacta- 
mente, saber calcular el valor de U a partir de cualesquiera 
valores fijos de Ry, Ri ++ 3 Can Come sed + 

La distribución probabilistica de las parámetros de cual- 
quier elemento se puede obtener experimentalmente anali- 
zando una gran partida de estos elementos. Dicha distri- 
bución resulta muy a menudo normel. Por eso, muchos 
investigadores parten de ello considerando, por ejemplo, 
que la resistencia del elemento representado en la fig. 
es una variable alaatoria normal p con la esperanza mate- 
mática Мр = 22 y cou За == 1,1 (recordemos que en un 
experimento es prácticamente imposible obtener un 
valor de p que difiera de Mp en más de 3c; véaso la fór- 
mula (20). 

El esquema del análisis es muy sencillo: para todo ele- 
mento se sortea primero el valor del parámetro y se calcula 
después, a partir de la formula (30), el valor de U. Este expe- 
rimento se repite N veces; obtenidos los valores Ch, Us, ... 
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-, Uy, Se toma aproximadamente 


N 
1 
MUR TÈ и, 
jam $ 


N N 
DU [2 Uy -— р (5 uj]. 
= pa 


Si М toma valores grandes, podemos sustituir el factor 
A/(N — 1) de la última fórmula por 1/N y entonces dicha 


Fig. 18 


férmula resulta un simple corolario de las férmulas (8) y (9). 
En la Estadistica Matematica se demuestra que para valores 
poqueños de N conviene emplear el factor 1/(N — 1). 

6.2. Ejemplos de análisis de la seguridad. Supongamos que 
queremos estimar el tiempo medio de duración de una pieza 
aceptando que conocemos el tiempo de duración de cada uno 
do sus elementos. 

Si consideramos que el tiempo 4; de duración de todo 
elemento es una magnitud fija, no habrá dificultad en calcu- 
lar el tiempo / de duración de la pieza. Por ejemplo, en la 
fig. 18 viene representada esquemáticamente una pieza que 
deja de funcionar cuando falla cualquiera de sus саа 
en este caso tenemos 


t= min (tas бө; feat t). (31) 


En cambio para la pieza con un elemento doblado repre- 
sentada esquemáticamente en la fig. 19, tenemos 


t = món Цан Las máx (у ta); to) (32) 


ya que si falla, por ejemplo, el elemento Ne 3, la pieza con- 
tinuará funcionando a base del elemento M 4, 
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El tiempo de duracién de cualquier elemento es, en reali- 
dad, una variable aleatoria ©,,,. Cuando decimos que el 
plazo de duración de una bombilla es de 1 000 horas, esta- 
mos hablando del valor medio МӨ de la variable Ө; todos 
comprendemos que una hombilla se fundirá antes mientras 
que otra (de la mísma partida) durará más. 

Si conocemos las densidades de distribución de Өш, 
para todos los elementos que componen la pieza, podemos 


Ш E 


Pig, 19 


calcular МӨ por el método de Montecarlo aplicando el proce- 
dimiento del punto 6.1. En efecto, para todo elemento pode- 
mos sortear el valor de la variable Өш); sea £a, este valor. 
Ahora, a partir de la fórmula correspondiente (31) o (32), 
podemos hallar el valor ¢ de la variable aleatoria O, Repi- 
tiendo este experimento un número (М) suficiente de veces, 
podemos aceplar que 


N 
H 
Me Ra ў ta 
del 
donde £, es el valor de £ obtenido en el f-ésimo experimento. 
Cabe subrayar que no es nada fácil determinar la distri- 
bución del plazo de duración Ө, de los elementos; por 
ejemplo, para elementos que funcionan durante largo tiempo 
se hace difícil la organización del experimento pues habrá 
que esperar a que falle una cantidad suficiente de estos 
elementos. 
6.3. Otras posibilidades del método. Los ejemplos expues- 
tos permiten ver que la metodología del análisis de la eali- 
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dad de las piezas se basa en une ideu muy simple, Es preciso 
conocer las características probubilísticas de todos los 
elementos que constituyen la pieza y saber calcular la va- 
riable que nos interesa como una función de los parámetros 
de estos elementos. La simulación permite entonces abarcar 
las propiedades aleatorias da los parámetros. 

Aparte de la esperanza matemática y de la varianza de 
la variable que nos interesa, la simulación permite obtener 
mucha más información útil. Supongamos, por ejemplo, 


2 о 1 15 24 32 7 19 7 1 2 
 ——= Sn — ee — — el 
a 2 3 4 5 Ф 7х 
Pig. 90 


que hemos obtenido una gran cantidad № de valores 1), 
Us ..., Un de la variable aleatoria U. Basándonos en 
estos valores, podemos encontrar aproximadamente la 
densidad de distribución de U. Esta cuestión concierne, en 
realidad, a la Estadística ya que se trata de la elaboración 
de resultados de experimentos (realizados, өз verdad, en 
una MCE). Por eso, nos limitaremos a considerar un ejemplo 
concreto. 

Supongamos que hemos obtenido un total de М = 120 
valores ІЛ, Us, . + ., Urgo de la variable aleatoria U com- 
prendidos todos en los límites 


1< Uj «65. 


Dividamos el intervalo 1 < z< 6,5 en 11 intervalos 
iguales de longitud Ax = 0,5 (о en otro número de interva- 
los no muy grande ni tampoco muy pequeño) y calculemos 
cuántos valores 1); corresponden а cada intervalo. Hemos 
representado estas cantidades en la fig. 20. 

Dividiendo entre № = 120 estas cantidades, obtendremos 
las frecuencias correspondientes a cada intervalo. En nuestro 
caso estas frecuencias son: 0,017, 0; 0,008; 0,12; 0,20; 0,27; 
0,14; 0,16; 0,06; 0,008 y 0,017. 

En cada intervalo construimos ahora un rectángulo de 
área igual а la frecuencia de U; correspondiente a dicho 
intervalo (fig. 21). En otras palabras, la altura de todo rectán- 


53 


gulo es igual a 1а frecuencia dividida entre Ах. La figura 
escalonada así obtenida se denomina Aistograma. 

El histograma ofrece una aproximación de Ia densidad 
desconocida de Ja variable aleatoria U. Por esto, el área del 
histograma comprendida entre, digamos, z = 2,5 y т = 5,9 
nos da un valor aproximado de la probabilidad 


P (2,5 < U < 5,5) = 0,95. 


El cálculo (experimento) realizado nos permite afirmar, 
por consiguiente, que con una prohahilidad aproximada- 
mente igual a 0,95 la variable 7 está comprendida en el 
intervalo 2,5 < U< 5,5. 

En la fig. 21 hemos representado, a título comparativo, 
la densidad de la variable aleatoria normal 2! con los pará- 
metros а = 3,85 y б = 0,88. Si partiendo de esta densidad 
calculamos la probabilidad de que 2) pertenezca al intervalo 
2,5<E « 5,5, obtendremos el valor 0,911) bastante próxi- 
mo al valor anterior. 


А) Expliquemos cómo se efectúan los cálculos, 
Ségun (14), tenemos 
4 EN ds. 
61/2 


2,5 


Р (2,5 «cU < 5,5) = 


Efectuemos en la integral la sustitución т — а = at. Tendremos 
entonces 


Yon 
P (25«ct 5,5) - ү. T at, 


Va 


1,54 y b= 


ü 


2,5—4 
= 1,88. 

Para calcular el valor de la última integral hay que recurrir а Ja tabla 

de la nsi llamada integral de probabilidades Ф (х) en la que vieoon 

las valores de la función 


donde n= 


А 
2 “ғ 
Фо--т- Ye ^at. 
| 


Asi encontramos 
P (2,5 « t! < 5,5} == 0,5 [0 (1,88) +0 (1,54)]— 0,94. 
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6,4. Observación. Lamentablemente son raros, por ahora, 
los casos en los que se aplican semejantes cálculos. Es difí- 
cil explicar esto. Posiblemente se deba a que los disefiadores 
ignoran esta posibilidad. 

Además, para poder aplicar este método al disefio de 
piezas es necesario estudiar previamente las características 


Y 2 3 4 5 6 x 


Fig. 21 


probabilisticas de los elementos que componen la pieza. 
Éste es un trabajo voluminoso. Es verdad, que, conociendo 
estas características, se puede estimar la calidad de cualquier 
pieza formada por dichos elementos. También se puede 
analizar cómo varía la calidad al sustituir un elemento por 
otro, 

Es de esperar que en los años próximos estos cálculos se 
hagan habituales. En cuanto a las características probabilís- 
ticas necesarias de los elementos, serán ofrecidas siempre por 
las empresas que los elaboren. 
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$7. 
ANÁLISIS DEL PASO 
DE NEUTRONES A TRAVÉS 
DE UNA PLACA 


Las leyes probabilísticas de la interacción entre una par- 
tícula elemental (neutrón, fotón, mesón, etc.) y la materia 
son conocidas. Como regla, es preciso conocer las caracterís- 
ticas macroscópicas —densidades, flujos, etc.— de procesos 
en los que participa una gran cantidad do estas partículas. 
Esta situación se asemeja a la que hemos analizado en el $ 5 
y en el $ 6 y resulta muy apropiada para la aplicación del 
método de Montecarlo. 

Posiblemente sea en la Física de neutrones donde se aplica 
con mayor frecuencia el método de Montecarlo. Analizaremos 
la variante más sencilla del problema sobre el paso de neu- 
trones a través de una placa. 

7.1. Planteamiento del problema. Supongamos que una 
placa homogénea infinita 0 < т < h es bombardeada por 
un flujo de neutrones de energía E, siendo el ángulo de inci- 
dencia igual a 90?. Al chocar con los átomos de la placa, los 
neutrones pueden esparcirse elásticamente o pueden ser 
absorbidos. Supongamos, para simplificar, que el neutrón 
conserva, al esparcirse, su energía y que cualquier cambio 
en la dirección del neutrón, provocado por el choqne con 
un átomo, es igualmente probable (esto último sucede en 
materias constituidas por átomos pesados). En la fig. 22 
hemos representado las distintas suertes que pueden correr 
log neutrones: el neutrón (a) atraviesa la placa, el neutrón 
(b) es absorbido y el neutrón (c) es reflejado. 

Es preciso calcular la probabilidad p* de que el neutrón 
atraviese Ja placa, la probabilidad p^ de que el neutrón 
sea reflejado y la probabilidad p° de que el neutrón resulte 
absorbido. 

a interacción de los neutrones con la materia se caracte- 
riza en este caso mediante dos constantes 2, y 2, que se 
denominan sección de absorción y sección de dispersión, ves- 
pectivamente. Los índices с y s provienen de las letras ini- 
ciales de las palabras inglesas «capture» (captura) y «scat- 
terings (dispersión). 
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La suma de estas secciones lleva el nombre de sección 
completa 


х= х, +5. 


Las secciones tienen el siguiente significado físico: al 
chocar el neutrón con un átomo de la materia, la probabilidad 


Lv 


Pig, 22 


de absorción es igual а 2/2 y la probabilidad de esparci- 
miento es igual а X,/2. 

El recorrido libre А del neutrón (o sea, su recorrido entre 
un choque y otro) es una variable aleatoria que puede tomar 
cualesquiera valores positivos con la densidad de proba- 
bilidad 

p (z) = ®е=?=, 


Es fácil ver que la densidad de la variable A coincide 
con la densidad (26) de la variable aleatoria т en el caso del 
flujo elemental de pedidos. Esta analogía con el punto 5.2 
nos permite escribir inmediatamente la fórmula para el 
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recorrido libre medio 


МА —1/Z 
y la fórmula para el sorteo de А 
1 
A= —.1n т. 


Resta explicar cémo se escoge la dirección aleatoria del 
neutrón cuando se produce el esparcimiento. Puesto que 
existe en el problema una simetría respecto al eje т, esta 
dirección se determina plenamente con sólo indicar el ángu- 
lo q entre la dirección de la velocidad del neutrón y el eje 
Oz. Se puede demostrar!) que la condición de que todas las 
direcciones son igualmente probables es en este caso equi- 
valente a la condición de que el coseno de este ángulo p = 
= соз q esté uniformemente distribuido en el intervalo 
(—1, 1). Tomando en la fórmula (24) a = —1 y b — 1, 
obtenemos la fórmula para el sorteo de н: 

н = 2y — 1. 

7.2. Esquema de análisis basado en la simulación de 
trayectorias reales. Supongamos que еї neutrón ha expe- 
rimentado el k-ésimo esparcimiento en el punto de abscisa zs 
correspondiente al interior de la placa y ha comenzado a 
moverse en la dirección рд. 

Sorteemos el recorrido libre 

М = —(t/Z) In y 


y calculemos la abscisa del punto de colisión siguiente 
(fig. 23) 
Жы = Tp БА. 

Comprobamos la condición de que la placa ha sido atrave- 

sada 
Tr > he 

Si esta condición se cumple, concluimos el análisis de la 
trayectoria del neutrón y agregamos uno al contador de las 
partículas que atraviesan la placa. En el caso contrario, 
comprobamos la condición de reflexión 


туы < б. 


1) La demostración viene en el punto 9.4. 
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Si esta condición se cumple, concluimos el análisis de la 
trayectoria y agregamos uno al contador de las partículas 
reflejadas. Si tampoco se cumple esta condición, o sea, si 
0< ть, SA, ello significará que el neutrón experimenta 


Vig. 32 


la (k + f)-ésima colisión en el interior de la placa; habrá 
que sortear de nuevo el «destino» del neutrón en esta colisión. 

De acuerdo con el punto 4.1 escogemos el valor siguiente 
de y y comprobamos la condición de absorción 


y< EE. 


Si esta última desigualdad se cumple, concluimos el ané- 
lisis de la trayectoria y agregamos uno al contador de par- 
tículas absorbidas. En el caso contrario, aceptamos que 
el neutrón experimente un esparcimiento en el punto de 
abscisa z,,,. Realizamos entonces el sorteo de la nueva 
dirección del neutrón 

н =29 4 


y volvemos a repetir todo el ciclo (tomando, claro está, 
otros valores para Y). 
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Hemos escrito todos los valores de ү sin indices sobren- 
tendiendo que cada valor de ү se emplea una vez solamente. 
Para analizar un eslabón de la trayectoria es preciso emplear 
tres valores de y. 

Los valores iniciales para cada trayectoria son 


zQ—0 y y, = 1. 


Después de haber analizado de esta forma N trayectorias, 
resultará que V* neutrones han atravesado la placa, que N- 
neutrones han sido reflejados y que N° neutrones han quedado 
absorbidos. 

Es evidente que las probabilidades buscadas son 
aproximadamente iguales a las razones 


м a N- мо 
Pea, Р я Жу pat, 


Еп la fig. 24 hemos ropresentado el esquema sinóptico 
del programa correspondiente a este análisis, El índice ј 
representa el número de la trayectoria y el índice Æ corres- 
ponde al número de la colisión (a lo largo de la trayectoria 
que se está analizando). 

Este método de análisis, a pesar de ser muy natural, 
tiene sus deficiencias. En particular, resulta difícil calcular 
por este método la probabilidad p* cuando ella es muy peque- 
ña. Pero precisamente con este caso nos encontramos al ana- 
lizar la protección contra Ja radiación. 

Existen modificaciones ingeniosas del método de Monte- 
carlo que pueden ser aplicadas también en estos casos. 

Detengámonos brevemente en la variante más sencilla 
del análisis que se basa en el empleo de los así llama- 
dos «pesos». 

7.3. Esquema de análisis basado en el empleo de pesos 
que sustituyen la absorción. Consideremos el mismo proble- 
ma sobre el paso de los neutrones. Supongamos que a lo 
largo de la trayectoria se desplaza un «tren» formado por un 
gran número w, de neutrones idénticos. Al producirse la 
colisión en el punto de abscisa z,, un promedio de we (DJE) 
de neutrones del «tren» será absorbido y un promedio de 
ша (Z,/X) de neutrones será esparcido. 

Agreguemos la cantidad wy (Z,/ È) al contador de particu- 
las absorbidas y continuemos el análisis del «tren» esparcido 
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Entrada de datos 


D 


{ 


k=0, граф; ө--4 


diua А? 


sí | ло 
иге] 
[1 
Му Мр 24 |е 


Er < 1/2? 


NS, NO 4-1 | sí | no 


Тау = ivt d T 


ES TAREA 
Ге] 


ai Tn | 


Cálculo de p*, p- y p? 


i 


Salida de resultados, fin 
del problema 


Fig, 26 


| Ng, = Ni Ll 
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aceptando que todo el «tren» restante se mueve en una misma 
dirección. 

Todas las fórmulas que hemos dado en el punto 7.2 se 
conservan. La ünica diferencia consiste en que, después 
de cada colisión, la cantidad de neutroues en el «tren» dis- 
minuirá, siendo 


шы = шу (Х,/2), 


ya que una parte del «tren», con un total de Wp (E¿/2) de 
neutrones, será absorbida; además, la trayectoria jamás 
podrá terminarse con la absorción. 

La variable wą suele denominarse peso de neutrón y, en 
lugar de un «tren» formado por ш, neutrones, suele hablarse 
de un neutrón de peso шд. El peso inicial wọ so toma, como 
regla, ignal а 1. Esto no contradice nuestra hipótesis de 
que el «tren» es grande, pues es fácil ver que wọ es un factor 
común de todos los pesos w, que se obtienen al analizar una 
trayectoria. 

El esquema sinóptico del programa correspondiente a 
este análisis viene represontado en la fig. 25. Como puede 
verse no resulta más complejo que el esquema de la fig. 24. 
Sin embargo, se puede demostrar!) que siempre es más con- 
veniente calcular p* por el método de este punto y no por 
el método del punto 7.2. 

7.4, Observación. Existe una gran variedad de proce- 
dimientos de cálculo que se basan en el empleo de distintos 
pesos. No podemos detenernos aquí en estas cuestiones. 

Señalemos solamente que el método de Montecarlo es 
aplicable a problemas mucho más complejos relacionados 
con partículas elementales. Por ejemplo, al caso en que el 
medio se compone de materias distintas y tiene una estructura 
geométrica arbitraria o al caso en que la energía de la partí- 
cula varía después de cada colisión. Así mismo, hace posible 
el análisis de otros muchos aspectos de los procesos nucleares: 
permite prever la posibilidad de que el átomo se desintegre al 
chocar con el noutrón y de que surjan neutrones nuevos; 
permite determinar el momento de inicio y la condición 
de mantenimiento de la reacción en cadena; etc. 


1) La demostración viene en el punto 9.5. 
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Entrado de datos 


| mr | 
П 
| | 
= -iny 


+ 
Thei ZR -H Agua 
Y 


| | шк = шд (24/58) | 


i D 


Гак = Һр Ni, Л ушр (2/2) | 
H y 
day NT рам =2y— d 
si no 7 
| mere Ld 
Cálculo de p*, p“ y p? 
t 
Salida de resultados, fin 
del problema 


Fig. 25 
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$8. 
CÁLCULO 
DE LA INTEGRAL DEFINIDA 


Los problemas considerados en los parágrafos 5, 6 y 7 
eran de contenido probabilistico por su esencia y la aplica- 
ción del método de Montecarlo en estos problemas parecía 
bastante natural. Ahora consideraremos un problema pura- 
mente matemático, el cálculo aproximado de la integral 
definida. 

Puesto que el cálculo de integrales definidas equivale al 
cálculo de áreas, podríamos emplear el método del punto 
1.2. Aquí explicaremos otro método, más efectivo, que 
permite construir diferentes modelos probabilísticos para 
resolver este problema por el método de Montecarlo y expli- 
caremos cómo escoger entre estos modelos los «más conve- 
nientes». 

8.1. Método de cálculo. Consideremos una función g (1) 
definida en el intervalo a < z< b. Es preciso hallar el 
valer aproximado de la integral 


b 
Te | g (2) dz. (33) 


Tomemos una densidad de distribucién рь (x) arbitraria 
definida en el intervalo (а, b) (o sea, una función p, (x) 
cualquiera que cumpla las condiciones (15) y (16)). 

Además de la variable aleatoria 5 definida con la den- 
sidad p; (т) en el intervalo (а, Б), necesitaremos emplear 
la variable aleatoria 


n= g (Spi (8). 
Según (18), tenemos 


b 
Mn = | te (/p:(2)1 pe(2)d2=1. 


Consideremos ahora № variables aleatorias т, Na, . - -. ma. 
independientes e idónticas, y apliquemos а la suma de estas 
variables el teorema central del límite del punto 2.4. La 
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férmula (21) se puede escribir en nuestro caso asi 
N — 
1 D: 
РПІ Жт—1}<з у EA} = 0,997. (34) 
jet 


Escojamos № valores &, Ej, .. ., Ёк. La última relación 
nos dice que siendo N suficientemente grande tendremos 


aga 
Y rd едім (35) 


Al mismo tiempo nos dico quo con una elevada probabilidad 
el error de la aproximación (35) no pasa de 
3V DwN. 

8.2. Cómo escoger el esquema de cáleulo. Como hemos 
visto para calcular la integral (33) se puede emplear una 
variable aleatoria E cualquiera definida en el intervalo (а, b). 
En todos los casos 


Mn = M le (Ep; ŒI = 4. 
Pero la varianza Dy y, por consiguiente, también la esti- 


mación (35) del error depende de qué variable $ tomemos. 
En efecto, 


b 
Юте М (19) — P | te Gym o) dz — P. 


Se puede demostrar! que esta expresión será mínima 
cuando pg (х) sea proporcional a | g (2) |. 

Por supuesto, no podemos escoger fórmulas muy com- 
plejas para ру (г) ya que entonces resultará laborioso sortear 
los valores de £. Pere es posible guiarse por esta recomenda- 
ción a la hora de escoger la función p; (2) (véase el ejem- 
plo del punto 8.3). 

En la práctica el método de Montecarlo no se aplica para 
el cálculo de integrales de tipo (33) ya que existen métodos 
más exactos basados en las fórmulas de cuadratura. Pero la 
Situación cambía si pasamos a considerar integrales múlti- 
ples: Jas fórmulas de cuadratura se hacen muy voluminosas 
mientras que ol método de Montecarlo no necesita casi 
modificaciones. 


1) La demostración viene en el punto 9.6, 
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8.3. Ejemplo numérico. Calculemos el valor aproximado de la 
integral 
m2 
i= | sen z dz. 
9 
El valor exacto de esta integral es 
n 
j pan z dz «е [ — cos 2144. 
ô 


Emplearemos para el cálculo dos distintas variables aleatorias Ё: 
una con la densidad constante p, (5) = 2/л (о sea, uniformemente 


distribuida en el intervalo (0, x/2)) y otra com la densidad lineal 
PE (2) = 82/n?. En la fig. 26 hemos representado ambas densidades así 


y у=8х 
т 
1 ye senx 
2 
YR 
1 ме 
o Wh x 
Fig. 26 


como la función integrando sen г. La figura permite ver que la den- 
sidad lineal está en mayor consonancia con la recomendación del pun- 
to 8,2 de quo es deseable que pg (x) y sen z sean proporcionales, Por 
ello, es de esperar que el segundo procedimiento dé un resultado mejor. 

(a) Soa pp (х) = 2/1 en el intervalo (0, 1/2). La iórmula del sorteo 
de E se puedo obtener de la fórmula (24) tomando a = 0 y b = 1/2: 


Es nyia 
La fórmula (85) da ” 
Im У) eng; 
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Sea № = 10. Tomemos para los valores de y las ternas de cifras de 
la tabla A (multipliséndolsa por 0,004). Los cálculos auxiliares apare- 
cen en la tabla 

Tabla 1 


TIE JE 


1 в 1 8 


Y 


o.oo sec cs od os [oss 0,612 |o a 


Es ЕТ 


ШЕТЛЕЛЕ 1,029 | 1,275 10,521 


0,247) ШЕСЕ 


sen ыт 


9, smo, 776 


El resultado de los cálculos es 
F 20,952. 
(b) Sea ahora py (т) = Вл/л?. Para sortear E empleamos la ecua- 
ción (23) 


(8z/nt) dz — y; 
9 
de aquí encontramos fácilmente que 


л 
E=3 Vr 

La fórmula (35) da А 
д? sent, 
BN к тен 

per 
Sea № = 10. Los números y se acogen igual que en el caso (a). 
Los cálculos auxiliares aparecen em la tabla 2. 


Ie 


Tabla 2 


sen 


$5 


lo, 680(0,936]0,908|0, 784)0 93710, 748. 0,863 | 0,754 | 0,698 | 0,868 


67 


El resultado de los cálculos es 
1 2 1,016. 


El segundo procedimiento, como lo habiamos previsto, da un re- 
sultado más exacto. 


8.4. Sobre la estimación del error. En el punto 8.2 ya 
hemos señalado que el error absoluto que se obtiene al calcu- 
jar la integral J mediante la fórmula (35) no puede práctica- 
mente sobrepasar la magnitud ЗУ Dn/N. Pero el error abso- 
luto real es, como regla, considerablemente monor. Por eso, 
para caracterizar el error suele emplearse en la práctica otra 
magnitud denominada orror probable 


Bor = 0,675Y DWN. 


El error ahsoluto real depende de los números aleatorios 
que se han empleado en el cálculo y puede suceder que sea 2 
6 3 veces mayor o varins veces menor que 5,,. Es decir, 8р. 
indica tan sólo el orden del error y no sus límites!). 


Volvamos al ejemplo del punto 8.3. Empleando los valores que 
aparecen en las tablas 1 y 2 y la fórmula para la varianza indicada en 
el punto 6.1, podemos calcular aproximadamente las varianzas Dy 
en ambos procedimientos ?). 

En la tabla 3 indicamos para ambos métodos los valore 
aproximados de las varianzas Dn los errores probables Spr calculados 


Tabla 3 
Método | On | Sor | Soa 
(a) 0,256 | 0,103 | 0,048 


0,016 | 0,027 


1) Sobre el error probable véase el punto 9,7. 
*) En el caso del método (a) tenemos 


10 10 
Dna = (5 (gen wag (E sni) |= 
{тт ~ 


LN 
=r (4,604 — 3,670) = 0,256, 
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a base de estos valores y los errores absolutos reales 6,4, que se han 


obtenido on los cálculos, 
Como vemos, бы es efectivamente del mismo orden que Spy. 


APÉNDICE 


#9. 
DEMOSTRACION 
DE ALGUNAS 
PROPOSICIONES 


En este parágrafo insertamos la demostración de algunas 
proposiciones enunciadas en los parágratos anteriores que o 
bien nos parecen un tanto voluminosas para el nivel que 
bemos mantenido өп la exposición o bien requieren cono- 
cimientos más profundos de la Teoría de las probabili- 
dades. 

9.1, Argumentación del método de Neumann para el 
sorteo de una variable aleatoria (punto 4.3). El punto alea- 
torio Г está uniformemente distribuido en el rectángulo abcd 
(fig. 27) cuyo área es igual a M, (b — а)). La probabilidad 
de que el punto Г aparezca por debajo de la curva y = 
= p (х) y no sea rechazado de la consideración, es igual a 
la razón entre las áreas 


b 
$ pla) аг 
Mo($—a)  My(b—a) 


Та probabilidad de que el punto Г aparezca por debajo 
de la curva у = р (т) y corresponda al intervalo a’ < т< 


En el caso del método (b) tenemos 
4 M SenEji2 4 10 sen Ey yz 
a 4 

бая yii [21 y) (2 2) ]- 


а 
= Fy (6.875 6,717) 0,016. 
3) Véase el punto 9.3, 
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< b’ es también igual a la razón entre las áreas 


Myg(b—a) ` 
Por consiguiente, entre todos los valores que hemos esco- 
gido para E, la parte de los valores que pertenecen al inter- 


YA 


ма = 


yep(x} 


4 


Pig. 27 


valo (a’, Ь') ез igual al cociente 


ы 
Í pi) dz pr 
Mb] Mba j p(z)dz 


> 
que es lo que se quería demostrar. 

9.2. Densidad de distribución de la variable $ = e + 
+ of (punto 4.4). Se supone que la variable 6 es una varia- 
ble normal con la esperanza matemática MC = 0 y 1а 
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varianza Di = 1 de modo que su densidad es 
жа 


п) Vut 


Para caleular la densidad de distribución de la variable 
Ç’, escojamos dos números arbitrarios ғ, <; х, y calcule- 
mos la probabilidad 


Рі « E <a} =P (z «a4 06 T) = 
=P (Spc). 


Por consiguiente, 
y ае a 
2 
Tu e * dz. 
Ys СҮЙ 
Realizando en la última integral la sustitución z' = a + oz, 
obtenemos 


Pla << 2) = 


Л t ЕКІСІ 
таз , 

== \ e ағ, 

Vino 2 

de donde resulta (véase (14)) que la variable 6” es una varia- 

ble normal сол los parámetros Mt’ = а у DỌ = o°. 
9.3. Distribución uniforme de los puntos (y’, ү”) еп 

el cuadrado (punto 4.5). Puesto que las coordenadas del 

punto (y', y”) son independientes, su densidad p (=, y) 

es igual al producto de las densidades 


р (ж, y) = py (2) ру (0). 


Cada una de estas densidades es igual idénticamente a 1. 
Por consiguiente, р(х, y) = 1 (para O«zz«t у 0< 
< у « 1). Pero esto significa precisamente que el punto 
(y', y”) está uniformemente distribuido en el cuadrado de 
dimensión 1. 

9.4. Selección de la dirección aleatoria (punto 7.1). La 
dirección se puede determinar mediante un vector de lon- 
gitud 1 que arranca del origen de coordenadas. Los extre- 
mos de estos vectores se encuentran en la superficie esférica 
de radio 1. Las palabras «cualquier dirección es igualmente 
probable» significan que el extremo del vector representa 


P(nctUxu 
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un punto aleatorio 2 uniformemente distribuido sobre la 
superficie esférica; es decir, la probabilidad de que © per- 


tenezca a cualquier elemento dS de la superficie es igual a = . 


xk 


+ 


Fig. 28 


Tomemos en la superficie esférica las coordenadas esféri- 
cas (ф, їр) con el eje polar Oz (fig. 28). Tendremos entonces 


dS == sen p dp dy (36) 


con la particularidad de que O < «p s x, 0 «y « 2n. 
Sea p (e, tp) la densidad del punto aleatorio (q, Y). 
Debemos tener 


р (Фф, 9) do d = dS/An; 
de aqui y de la relación (36) se deduce que 


sen p 


Al т 5 (37) 
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A partir de la densidad conjunta de ф у xp, podemos calcular 
fácilmente tas densidades de estas variables: 


2л 

рае) = (рер de St, (38) 
0 

рь) = | pir. dou. (38) 


9 


La igualdad р (q, Y) = Po (9) Po (4) demuestra que p y 
y son independientes. 

Es obvio que y está uniformemente distribuida en el 
intervalo (0, 2л) y, por eso, la fórmula para el sorteo de ф es 


Y = 2ny. (40) 
La fórmula para el sorteo de p se obtiene de la ecua- 
ción (23): 


1 
- | зеп хіт= 7, 


оча 


de donde resulta 
cos q = 1 — 2%. 


Las fórmulas (40) y (41) permiten escoger (sortear) la di- 
rección aleatoria. Por supuesto, los valores de y en estas 
fórmulas deben ser distintos. 

La ültima fórmula del punto 7.4 se obtiene de lo fórmula 
(41) si se toma 1 — y en lugar de y; pero estas variables tienen 
la misma distribución. 

9.5. Ventaja del método de análisis basado en el empleo 
de pesos (punto 7.3). Consideremos las variables aleatorias 
v y v' iguales a la cantidad (peso) de neutrones de una misma 
trayectoria que atraviesan la placa, con la particularidad 
de que la primera variable corresponde al método del punto 
7.2 y la segunda, al método del punto 7.3. 

Según el contenido físico del problema tenemos 


Mv = Mv’ = р". 


(Véase Ja demostración rigurosa de esta alirmación en el 
libre [3].) 
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Puesto que v puede tomar sólo dos valores, O y 1, la dis- 
tribución de v es 
1 0 
Ves eee 


Tomando en cuenta que v? = v, os fácil ver que Dv = 
= р" — (p9* Es evidente que la variable y” puede tomar, 


aparte del valor 0, un conjunto infinito de valores wy = 1, 


Wa. Wer -s Up, ... Por eso, su distribución es 
vin (то 
do 9 Yas e+ 4» -+q 


Los valores concretos de g,, como veremos, no nos harán falta. 
Cualesquiera que sean podemos representar la varianza en 
la forma 


DY = д wga — (РУ. 


Observando que todos los w, < 1 y que Djwpg = Mv’ = 
' c=) 


= p*, obtenemos la desigualdad Dv’ < p* — (p*)? = Dv. 

El hecho de que la varianza de v' sea siempre menor que 
la varianza de v significa que, al calcular p*, el método del 
punto 7.3 dará siempre mejores resultados que el método 
del punto 7.2. La misma conclusión podemos hacer en cuanto 
al cálculo de p- y también al cáclulo de p° si la absorción 
no es muy grande. 


9.6. Selección óptima de E (punto 8.2). En el punto 8.2 
hemos obtenido la expresión para la varianza Dy. Para 
hallar el mínimo de esta expresión, siendo arbitraria la selec- 
ción de p; (х), recurriremos a la siguiente desigualdad 
conocida del Análisis 


ri lu (z) to! afc} w(x) dz- | 0 (2) dz. 


Tomando и = g (zVV P; (2) y v = V Pela), obtenemos de 
esta desigualdad 


b 
Пеене) E 
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Por consiguiente, 
b 
т>Г еда f- r. аз) 


Reste demostrar que la cota inferior se alcanza cuando 
pa (2) es proporcional a | к (x) |. 
ea 


pz) = LLL... (43) 
AGILA 


а 
Es fácil calcular que para esta densidad р; (z) se tiene 


b b 

М » q2 

j tevin @lae=[ | Le (2) 1а] 

a & 
de modo que la varianza Dn resulta efectivamente igual al 
segundo miembro de (42). 

Observemos que es prácticamente imposible emplear en el 
cálculo la densidad «más conveniente» (43): deberíamos 
b 


conocer para ello el valor de la integral j lg (2) | dz; pero el 


a 
cálculo de esta integral es, de hecho, un problema equiva- 
lente al probloma que tenemos planteado acerca del cálculo de 


la integral | g (x) dz. Por eso, nos hemos limitado a dar la 


2 
recomendación indicada en el punto 8.2. 

9.7. Definición del error probable (punto 8.4). Sea £ la 
variable aleatoria normal definida en el punto 2.3. Es fácil 
calcular que para г = 0,6750 y cualesquiera a y 0 se tiene 


афт 
{ Py (2) dz = 0,5. 
a-r 

De aquí so deduce que 
P {&—а|<г}=рР (18—a]>r) = 0,5, 

о sea, que las desviaciones шауотез чие г у las desviaciones 

menores que ғ son igualmente probables. Por eso, la mag- 

nitud г se donemina error probable de la variable aleatoria E. 
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En el punto 8.1 hemos realizado-el cáleulo aproximado de 
la variable normal p = (1/N) (n, + n, +... + n, Re- 
sultó que su esperanza matemática es а = Mp = / y que 
su varianza es 0? = Dp = Dn/N. Por eso, el error probable 
de la variable р es aproximadamente igual a 0,675V Dn/N. 


8 10. 
SOBRE LOS NÜMEROS 
SEUDOALEATORIOS 


Los algoritmos que se emplean para obtener números 
seudoaleatorios tienem, en su mayoría, la forma 


dau Е (ч). (44) 
Si se ha dado el número inicial yy, todos los números siguien- 
tes ү, үз... se calculan aplicando la fórmula (44) 
con k — 0, 1, 2, ... Ambos algoritmos considerados en el 


punto 3.3 son también de la forma (44) sólo, en lugar de 
señalar la forma analítica de la función y = F (xj, se ha 
indicado el conjunto de operaciones que debo efectuarse con 
el argumento z para obtener el valor de y. 

10.1. Cómo debe escogerse la función F (ж). El ejemplo 
que sigue explica еп qué consiste una de las dificultades 
principales que se plantea al escoger F (т). 


Ejemplo. Demostremos que la función y = /'(z) cnyo gráfico 
aparece en Ja fig. 28 no se puede emplear para obtener, a partir de la 
formula (44), números seudoaleatorios, 

En efecto, consideremos en el cuadrado 0 < x < 1, 0 «y «1 
los puntos de coordenadas 

s т)» (Yas yas Qs POs s 
donde y; = F (ру), ya = F (Ys), Ya = F (үз), ... Todos estos puntos 
$e encueutran sobre la curva y = F (x) y esto, por supuesto, está muy 
mal ya que los puntos aleatorios verdaderos deben embrir uniiorme- 
mente todo el cuadrado (véase el puuto 9.3). 


Este ejemplo pormite ver que el empleo de la función 
у = F(z) on la fórmula (44) puede resultar exitoso sólo 
Si el gráfico de esta función recubre con bastante densidad 
todo el cuadrado. 

Este es el caso, por ejemplo, de la función 


y = {ez}, (45) 
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donde g es un nümero muy grande у {z} representa la parte 
fraccionaria del nümero z. En la fig. 30 aparece el gráfico 


Fig. 20 


de la función (45) correspondiente a g — 20. El lector podrá 
hacerse una idea de lo que representará este gráfico si g = 517. 

10.2. Método de congruencias (método de residuos). El 
algoritmo que más se emplea para obtener nümeros seudo- 
aleatorios fue propuesto por D. Lehmer. Se basa en la fun- 
ción (45), pero suele construirso de una manera algo distinta 
que permite realizarlo fácilmente en las MCE. 

Se determina una sucesión de nümeros enteros m, toman- 
do m, = 1 y calculando todos los números sucesivos ту, 


Ma, ... а partir de una misma fórmula 
ткы = 5m, (mod 2%), (46) 
donde k = 0, 1, 2, ... Basándose en los números Me, se 


calculan los números seudoaleatorios 
Y = 2 ту. (47) 
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La fórmula (46) significa que el número m,+, es igual al 
resto que se obtiene al dividir 517m, entre 2%. 120 la teoría 
de las congruencias (véase cualquier texto de Teoría de los 


y ={20x} 


Fig. 30 


nümeros) este resto lleva el nombre de residuo positivo mó- 
dulo 2% mínimo. Por eso, el algoritmo suele denominarse 
método de congruencias o método de residuos. 

Es muy fácil realizar las fórmulas (46) y (47) en las 
MCE que funcionan con nümeros de 40 órdenes aplicando 
la operación de multiplicación con el nümero duplicado 
de órdenes!) hay que emplear las cifras correspondientes а 
los órdenes menores del producto. 

El período de la sucesión m4 es igual a 2%, 

En el caso de las MCE que funcionan con numero de 36 órdones, se 


toman en las fórmulas (46) y (47) los números 5% у 2% en lugar de 517 
y 240; el período de la sucesión ту es igual entonces a 2%, 


1) La MCE Strelá (véase el punto 3.3) no puedo 
realizar esta operación. 


yo 


Tablas 


Tabla A. 400 cifras aleatorias 1) 


1,0423 
1,8818 


Las cifras aleatorías imitan los valores de la variable aleatoria con la distribución (22) {véase e] punto 3,1). 
Los valores normales imitan los valores de la variable aleatoría normal (gaussiana) E con los parámetros a = 0 
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